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I review results recently obtained within the Hamiltonian approach to Yang-Mills theory in 
Coulomb gauge. In particular, I will present results for the ghost and gluon propagators and 
compare these with recent lattice data. Furthermore, I will give an interpretation of the inverse 
of the ghost form factor as the dielectric function of the Yang-Mills vacuum. Our ansatz for 
the vacuum wave functional will be checked by means of functional renormalization group flow 
equations, which are solved for the gluon energy and the ghost form factor. Finally, we calculate 
the Wilson loop for the vacuum wave functional obtained from the variational approach, using a 
Dyson equation. 
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Hamiltonian approach to Yang-Mills theory in Coulomb gauge 



1. Introduction 



The usual language of quantum field theory is the functional integral framework. However, from 
ordinary quantum mechanics we know that for many non-perturbative studies the solution of the 
Schrodinger equation is much simpler than calculating the corresponding functional integral. Al- 
though in quantum field theory the regularisation and renormalisation procedures are much better 
understood in the functional integral formalism, for non-perturbative investigations like the infrared 
sector of Yang-Mills theory, the canonically quantised operator formalism seems to have certain 
advantages over more traditional field theoretical approaches if it comes to the computation of 
physical observables. In recent years there have been many activities in studying the infrared sec- 
tor of QCD in Coulomb gauge. The use of the Coulomb gauge is advantageous since this gauge is a 
so-called "physical gauge": The gauge degrees of freedom can be directly removed and Gauss' law 
can be explicitly resolved. For example, the confining potential between static colour sources can 
be extracted much more easily than e.g. in Landau gauge. In this talk I report on results obtained 
in recent years by our group by a variational solution of the Yang-Mills Schrodinger equation in 
Coulomb gauge. The plan of my talk is as follows: 

After these introductory remarks I will briefly review the basic ingredients of the Hamiltonian 
approach to Yang-Mills theory in Coulomb gauge. In Sect. 3 I will present some results obtained 
from a variational solution of the Yang-Mills Schrodinger equation in Coulomb gauge and compare 
them with recent lattice data. There I will give a physical interpretation of the ghost form factor in 
Coulomb gauge. Our variational ansatz for the Yang-Mills wave functional will be checked in Sect. 
4 by means of the functional renormalisation group flow equations, which will be solved assuming 
ghost dominance. In Sect. 5 and 6 I will present some applications of our approach. Thereby I 
will focus on the calculation of the topological susceptibility and of the Wilson loop by means of a 
recently proposed Dyson equation. 



2. Hamiltonian approach to Yang-Mills theory in Coulomb gauge 



Standard canonical quantisation of Yang-Mills theory in Weyl gauge Aq = leads to the Hamilto- 
nian 

H= l -J(Yl 2 +B 2 ), (2.1) 

where Yl = 8/ id A is the momentum operator and B" the non-Abelian magnetic field. Due to the 
Weyl gauge Gauss' law escapes from the equation of motion and has to be implemented as a 
constraint on the wave functional 

DI%) = P|V) ■ (2-2) 

Here D = d +gA,A ah = f acb A c is the covariant derivative in the adjoint representation of the gauge 
group and p denotes the colour charge of the matter fields. Implementing the Coulomb gauge dA = 



Gauss' law ( |2.2| ) can be explicitly resolved resulting in the Coulomb gauge fixed Hamiltonian 

H = \J [r l n L m L +B 2 ^ + 1- j V 1 (p+p dyn ) JF{ P + Pdyn ) , (2.3) 
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where H 1 - denotes the transversal momentum operator, J = Det(-Dd) is the Faddeev-Popov de- 
terminant and p a dyn = —Af'Yl h k L is the colour charge of the gauge field. Furthermore 

F ab (x,y) = (a,x\(-Dd)-\-d 2 )(-Dd)- l \b,y) (2.4) 

is the so-called Coulomb kernel. Its vacuum expectation value (F) defines the non-Abelian Coulomb 
potential. In the gauge fixed theory the matrix elements of an operator f?[A,IT] between states of 
the Yang-Mills Hilbert space are defined by 

(Yi\0[A,U}\Y2) = f DA x J[A x ]\fi[A}0[A,n]\ir 2 [A} , (2.5) 

where the restriction to the integration over the transversal gauge fields and the Faddeev-Popov 
determinant in the integration measure arise from the implementation of the Coulomb gauge. 
In Ref. ^ a variational solution of the Yang-Mills Schrodinger equation has been accomplished 
using the following ansatz for the vacuum wave functional 1 



VA[A]=7[A]- 1 / 2 exp ? 
where co(x,y) is a variational kernel, which is determined by minimising the energy 



-r / AcoA 



(2.6) 



(y\y) 



For the wave functional ( |2.6| ) the static gluon propagator is given by 

(AA) = (2fl)) _1 (2.8) 

implying that co represents the gluon energy. In Ref. [|I|] the energy (H) was calculated up to two 
loops. Furthermore the gap equation 8(H) / 8(0 = was converted into a set of Dyson-Schwinger 
equations. 

3. Results 

An infrared analysis, Ref. [f|], of the Dyson-Schwinger equations shows that the infrared exponents 
of the gluon and ghost propagators defined by 

(0{k)~k- a , d(k)=k 2 ((-Ddy l ) ~k-P (3.1) 

satisfy the sum rule 

a = 2/3 - 1 (3.2) 

and allow for two solutions 

i) p -0.796(0.85) ii) = 1.0(0.99) (3.3) 



'in Refs. |g] a pure Gaussian ansatz was used. Furthermore, in Ref. [g| the Faddeev-Popov determinant was 
completely ignored, while in Ref. [^] it was included in the kinetic part of the Yang-Mills Hamiltonian only. For a more 
general discussion of the ansatze for the vacuum wave functional, see Ref. [M 
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Figure 1: Ghost form factor (left panel) and gluon energy co(k) (right panel). 



which were also produced by the numerical solutions obtained in Refs. [[I]] and respectively. 
The corresponding numerically obtained infrared exponents are given in the brackets. In the nu- 
merical solution the horizon condition 

d-\k = 0)=0 (3.4) 

was explicitly built in. 2 Figures |] and || show the ghost form factor and gluon energy as functions 
of the momentum for the solution ii). At large distances the gluon energy rises linearly with the 
momentum as expected from asymptotic freedom, while it is infrared divergent at small momenta. 
Solution ii) gives rise to a strictly linearly rising static colour potential shown in fig. ^[ The running 
coupling constant extracted from the ghost-gluon vertex obtained for the solution ii) is shown in 
fig. H (left). It is the solution ii) shown in figs, [j] and ||, which is also supported by the lattice 
data obtained in Ref. [JSp _ Remarkably, the lattice gluon energy can be perfectly fitted by Gribov's 
formula 

with the energy scale M « 0.8 GeV. The lattice calculations carried out in Ref. ^ differ from 
previous lattice calculations in two respects: the residual gauge invariance left after implementing 
the Coulomb gauge has been explicitly fixed and the scaling violations have been eliminated giving 
rise to a strictly multiplicatively renormalisable static gluon propagator. For more details we refer 
to Ref. [§. 

In Ref. [Q] it was shown that the inverse of the Coulomb gauge ghost form factor d~ l (k) can be 
interpreted as the dielectric function of the Yang-Mills vacuum 

e(jfe) =d~\k) , (3.6) 



which is shown in fig. |3| (right). By the horizon condition (3^) the dielectric constant vanishes in 



the infrared e(k = 0) = 0, implying that the Yang-Mills vacuum is a perfect dia-electric medium, 
in which by Gauss' law V(eE) = p free colour charges cannot exist and thus have to be confined 
in colourless states. This is precisely the picture underlying the MIT bag model, which assumes 
that £ = in the vacuum while e = 1 inside the hadronic bags. The magnetic analogue of a perfect 



2 In D = 3 + 1 there are also so-called subcritical solutions H7||, satisfying d 1 (0) ^ 0. 
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Figure 2: Static gluon propagator 1 / (2<x>(k)) compared to the lattice results (left panel) and static quark 
potential (right panel). 
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Figure 3: Running coupling constant (left panel) and dielectric function of the Yang-Mills vacuum (right 
panel). 

dia-electric medium is a superconductor for which the magnetic permeability p: defined by B = \iH 
vanishes (B-magnetic field, //-induction). In the usual notion of duality, meaning the interchange 
between electric and magnetic fields and charges, a medium with e = is a dual superconductor. 
Therefore, the Gribov-Zwanziger confinement scenario assuming d~ 1 (k = 0) =0 implies by the 
identity < J3-6| ) the dual Meissner effect in the infrared. 

4. Functional renormalisation group flows 

Our ansatz for the Yang-Mills vacuum wave functional can be tested by comparison with results 
from functional renormalisation group flows (FRG). The basic idea of the hamiltonian FRG is to 
add an infrared cut-off term quadratically in the quantum field <p 

AS k l$} = ^jd 3 p$(p)R k (p)$(-p) (4.1) 

to the Euclidean action, which cuts off momentum modes of the field with momenta p < k, but 
leaves the theory unchanged for momenta p > k. At large cut-off scales k — ► Aj/y the theory is well 
under control due to asymptotic freedom and perturbation theory can be applied. In turn, for k — ► 
one recovers the full theory. The flow of the theory with the cut-off scale k is described by the 
renormalisation group flow equation for the infrared regularised effective action r\[0]. For cut-off 
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Figure 4: FRG-hows for the propagators. 
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Figure 5: FRG-flows for the propagators assuming ghost dominance and dropping the tadpoles. 



terms (4.1) the flow equation reads 

1 



d,r*[0] 



Tr 



<2) 



[<!>]+ R k ) d t R k 



d t = k 



dk 



(4.2) 



(2) 

where F K k } [0] + Rk is the inverse propagator of with 



8 2 r k [<i>] 

8<p8<p 



(4.3) 



The flow equation ( |4.2| ) entails the evolution of the IR-regularised effective action from k — > < 
where r^[0] coincides with the bare action S[0] to /c — » where r\[0] is the full effective action. 



Here we only are interested in the flow of the propagators, which is obtained from (4.2) by taking 
the second functional derivative with respect to the fields 

8 2 



dr {2) 



5050 2 



-Tr 



if'M+Rk) d t R k 



-l 



(4.4) 



0=0 



This equation is diagrammatically illustrated in fig. f| for the propagators of Yang-Mills theory 
and is structurally similar to a Dyson-Schwinger equation except that the infrared regulator d t R k 
enters the loops and all vertices and propagators are fully dressed. For the Hamiltonian flow of 
Yang-Mills theory in Coulomb gauge the fields involved are the transversal gauge field and the 



ghost field. Thus the right-hand side of the flow equation ( |4.4[ ) receives contributions from ghost 
and gluon loops. We assume ghost dominance and keep only the contributions from the ghost loop 
to the right-hand side of the flow equation (fi~4|). The resulting flow equations for the gluon and 
ghost propagators are diagrammatically illustrated in fig. |[ With our choice of the vacuum wave 
functional the action is given by 
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Figure 6: FRG-flow of the ghost form factor d(p) (left panel). Topological susceptibility obtained in the 



variational approach from solution ii) ( p.3| ) as function of G c / a (right panel). The two horizontal lines limit 
the range of the lattice results. 



S[A,c,c] = J ACOA + J c(-Dd)c . (4.5) 

We solve the FRG flows shown in fig. || for the gluon energy (o(p) and the ghost form factor d(p) 
using the following regulators 

( k 2 p 2 \ ( k 2 p 2 \ 

gluon: R k (p) ~ p exp f - j , ghost: R k (p) ~ p 2 exp ( -j - \ , (4.6) 

where we have suppressed the tensor structure, and perturbative initial conditions at large cut-off 
scale k = A 

(Oa{p) = p + const , d\{p) = const (4.7) 

down to some minimal momentum cut-off scale k m i n . Fig. ^ illustrates the renormalisation group 
flow of the ghost form factor. As the cut-off scale k is reduced the ghost form factor as function 
of the momentum p gets infrared enhanced and eventually becomes infrared divergent as the cut- 
off scale k becomes very small, i.e. the horizon condition emerges when the infrared cut-off is 
removed. This is nicely seen in fig. [7] where the ghost form factor is shown at the minimum cut-off 
k m i„ as a function of the momentum p. It is also seen that the infrared exponent, i.e. the slope 
of the curve dk min {p), does not change as the minimal cut-off scale is lowered. Fig. || shows the 
corresponding result of the integration of the flow equation for the gluon energy (Ok mjn {p)- The 
infrared exponents of the gluon energy and the ghost form factor satisfy the sum rule ( |3.2[ ) found 
from the Dyson-Schwinger equation of the variational approach. It can be shown that the FRG 
results coincide analytically with those of the variational approach either for optimised regulators 



[|10p, or if the ghost tadpoles are included. From analogous studies in the Landau gauge [ |11[ ] we 
expect that for general regulators the infrared exponents are smaller. Indeed this is the case for the 
approximation illustrated in fig. |5] and with the regulators (4.6). This can be seen from figs. [7] and 



|8[ The gaps between the solutions in figs. [7] and |8] and the difference in the exponents gives an 
estimate for the systematic error of the present approximation. A detailed analysis will be provided 
elsewhere. 



7 



Hamiltonian approach to Yang-Mills theory in Coulomb gauge 



10000 




1000 



le-08 le-06 0.0001 0.01 1 100 
p [GeV] 




p [units DSE] 



Figure 7: Ghost form factor df,{p) (left panel). Forsake of comparison the result of the variational approach 
is also shown (right panel). 
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Figure 8: Gluon energy (£>k{p) obtained from the solution of the RG- flow equation (left panel). For sake of 
comparison the result of the variational approach is also shown (right panel). 



5. Physical applications 

An important quantity for hadron physics is the topological susceptibility % defined by 

g 2 



X 



d 4 x(0\q(x)q(0)\0) 



(5.1) 



where q(x) is the topological charge density. It reflects the anomalous U(l) symmetry breaking in 
QCD and can be extracted from the 7]' mass through the Witten-Veneziano formula 



m*, + m 2 



F n = 93 MeV. 



(5.2) 



This quantity vanishes to all orders in perturbation theory and is therefore perfectly suited to test 
non-perturbative methods. % is a manifestation of the 6 —vacuum and can be easily evaluated in 
the Hamiltonian approach. Adding the topological term Jtf t0 p = Q J d 4 x q{x) to the Yang-Mills 
Lagrangian shifts the momentum by IT — > IT — d f-j B and the Hamiltonian of the 6— vacuum reads 



"°=2 



(53) 



from which one derives the following expression for the topological susceptibility (V is the spatial 
volume) 
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1 dd 2 



(5.4) 

=o 



Using the Yang-Mills wave functional obtained from the variational principle for the solution ii) 
(3.3) and restricting the intermediate states in eq. (5.4) to 3-gluon states (on top of the vacuum) one 
finds [12] the results shown in fig. ||| (right panel), where % is plotted as function of the ratio o c j o. 
a c is the Coulomb string tension, which is extracted from the static colour potential and used to fix 
the scale in the variational calculations while a is the Wilsonian string tension extracted from the 
Wilson loop. Lattice calculations indicate that this ratio is in the range of o c /o = 1.5. For such 
ratios the value obtained for % is somewhat larger than the quenched lattice results. 
A crucial test for a wave functional would be the calculation of the Wilson loop, the order parame- 
ter of the Yang-Mills theory. The spatial Wilson loop can in principle be directly calculated in the 
Hamiltonian approach, once the vacuum wave functional is known. However, the practical calcu- 
lation is rendered difficult by the path ordering. A quantity more easily accessible is the 't Hooft 
loop, the disorder parameter of Yang-Mills theory. In Ref. JT3| ] the 't Hooft loop was calculated 
using the wave functional corresponding to the solution ii) ( p.3[ ) and a perimeter law was found, 
which is characteristic of the confinement phase. For further details we refer the reader to Ref. 

S 



6. Wilson loop from a Dyson equation 

Although the temporal Wilson loop is dual to the spatial t'Hooft loop and a perimeter law in the 
latter implies an area law in the former, one would appreciate an explicit calculation of the Wilson 
loop and observe the emergence of the area law. Recently, a Dyson equation has been proposed 
in the context of supersymmetric Yang-Mills theory, which, at least in an approximate fashion, 
takes care of the path ordering Jl4]]. It has been applied to the temporal Wilson loop in non- 



supersymmetric Yang-Mills theory [15] 



This Dyson equation applies to trapezoidal loops W(S,T;L) and basically sums up the ladder 
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(a) 




Figure 10: (a) Processes involved and (b ) processes ignored in the Dyson equation ( 6.1). 



diagrams shown in rig. g(a), i.e. the gluon exchange between one pair of opposite paths. This 
Dyson equation is diagrammatically shown in fig. ^(b) and is analytically given by 

5 T 

W{S,T-L) = l+g 2 C 2 jds J dtD^(x(s)-s(t)) 2 ^W(s,t;L) . (6.1) 



This equation suffers from the following limitations [|l6|] : 

1 . Since we include only one pair of paths (see fig. |l(](a)) the Dyson equation ( |Q| ) is restricted 
to strongly asymmetric loops. However, it is irrelevant whether the loop is a temporal or 
spatial one. 

2. When one of the two parallel temporal lines is set to zero, the trapezoidal loop degenerates 
to a triangle-shaped loop. In this case the Dyson equation (FTTh yields for the Wilson loop 



W (S,T =0;L) = 1 (6.2) 

which, in general, is certainly incorrect. This wrong boundary condition is not surprising, 
since setting T = or S = contradicts the assumption S, T S> L in the Dyson equation (|Q|). 



3. The equation is bounded to produce strict Casimir scaling, which is known to occur only in 
the intermediate distance regime. 

4. The Wilson loop is a gauge-invariant object. However, the right-hand side of the Dyson 
equation ( |6T| ) depends via the gluon propagator on the gauge. In fact, the gluon propagator 
can be defined only in a gauge-fixed theory and vanishes if the gauge is unfixed. 

5. Finally, the Wilson loop is renormalisation group invariant, while the right-hand side of the 



Dyson equation (6.1) is not, except for the temporal Wilson loops in Coulomb gauge (which 



we will consider at first). 
In Coulomb gauge the temporal gluon propagator has the form 

g 2 D$ ) (x,y) = -8 ah V c (\x-y\)8(x -y°)+P ah (x,y) (6.3) 

where Vc(l x — y|) is the so-called non-Abelian Coulomb potential. At small distances it has the 
ordinary Coulombic ~ l/r behaviour, while it rises linearly at large distances ~ a c r, with a c 
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L[fm] L [fm] 

Figure 11: Left panel: The Wilsonian potential V(L) obtained from the static gluon propagator < ft3i ) and 
the perturbative potential V pert (L) . Right panel: The full potential minus its perturbative part. 



being the Coulomb string tension, which is larger than the Wilsonian string tension o. The non- 
instantaneous part P(x,y) is assumed to lower o c towards o. If we ignore the screening part P(x,y), 



the gluon propagator is instantaneous and the Dyson equation ( p.lD applies only to rectangular 
loops W(T;L) =W(S = T,T;L) 



W(T;L) = 1 -C 2 V C {L) J dtW(t;L) 
b 



(6.4) 



This equation contains the correct boundary condition W (T = 0;L) = 1 and can be solved analyti- 
cally, yielding 

W(T;L) = exp (-C 2 V C (L)T) . (6.5) 

Within the approximation used for the gluon propagator, we have correctly obtained an area law. 
It is clear why in this case the equation produces the correct result: The processes neglected in the 
Dyson equation (6.1) (see fig. (|T0|(b)) do not exist for an instantaneous gluon propagator. 



For arbitrary (non-instantaneous) gluon propagators, the Dyson equation (6.1) can be converted 
into a one-dimensional Schrodinger equation 



I 2 
dr 



a 2 

9n{r) = —£<Pn(r) 



with the variable r = and the Schrodinger potential given by 

U(r) = -g 2 C 2 L 2 D(L 2 (l+r 2 )) . 

The Wilsonian potential can be obtained from the "ground state energy" 

Oo(L) 



V(L) 



lim -lnW(T,T;L) 



+ const 



(6.6) 



(6.7) 



(6.8) 



Applying the Dyson equation to the spatial Wilson loop in Coulomb gauge and using the static 



transversal gluon propagator (3.5) augmented by its anomalous dimension (derived in Ref. [jl 
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to make the gluon propagator well-defined in coordinate space, one finds the Wilsonian potential 
shown in fig. [flja). When one subtracts from this potential the perturbative one, a linearly rising 
potential is left, implying an area law for the spatial Wilson loop, see fig. [TT|(b). 
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